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Abstract 

Radiative corrections to the decay rate of n = 2 states of hydrogenic ions are calculated. The 
transitions considered are the Ml decay of the 2s state to the ground state and the E1(M2) decays 
of the 2p]/2 and 2j> 3 / 2 states to the ground state. The radiative corrections start in order a(Za) 2 , 
but the method used sums all orders of Za. The leading a(Za) 2 correction for the El decays 
is calculated and compared with the exact result. The extension of the calculational method to 
parity nonconserving transitions in neutral atoms is discussed. 
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I. INTRODUCTION 



Radiative corrections to decay rates in atoms and ions have not been as thoroughly 
studied as other kinds of radiative corrections, such as the Lamb shift Q and corrections to 
hyperfme splitting (hfs) . An exception is the case of the exotic atom positronium, where 
differences between the lowest-order decay rates and experiment of -0.6 and -2.2 percent are 
present for parapositronium j^J and orthopositronium respectively. In both cases the 
bulk of the difference is accounted for by one- loop ^ivec—S 3, which mter 
in order a with large coefficients, and the present theoretical interest has advanced to the 
level of two-loop radiative corrections [8,0]. 

For other hydrogenlike atoms, theoretical work on one-loop corrections to Ml decays has 
been carried out in Ref. 111. These papers established that, unlike positronium, the 
order a correction has a vanishing coefficient, but did not calculate the actual correction, 
which enters in order a(Za) 2 . A calculation for El decays of the a(Za) 2 \n(Za) correction 
has been carried out in Ref. ITJ, and is in disagreement with another calculation associated 

n 

with the experimental determination of the Lamb shift |13( ■ This situation will be discussed 
further in the conclusion. 

It is the purpose of the present paper to calculate radiative corrections for the hydrogen 
isoelectronic sequence using methods that treat the electron propagator exactly. In addition, 
a perturbative calculation for El decays through order a(Za) 2 is carried out and compared 
to the exact result. While of intrinsic interest, development of these techniques should also 
aid in the evaluation of radiative corrections to parity nonconserving transitions in atoms, 
as will be discussed in the conclusion. 



II. LOWEST-ORDER CALCULATION 



While the first calculations of the decay rate of hydrogen date back to the beginning of 
quantum mechanics, fully relativistic calculations needed for calculations of highly-charged 

n 

hydrogenlike ions were first carried out in the early 1970's [14]. We briefly present the theory 
here using techniques that will be extended to the radiative correction case. We want to use 
the fact that a decay rate can be related to the imaginary part of the energy through 

T = -2Q(E). (1) 
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This is the approach taken by Barbieri and Sucher The one-photon decay rate is 

connected through this formula with the self-energy of an electron in a state v, which will 
be chosen here to be 2si/ 2 , 2py 2 , or 2p 3 / 2 . We define this self-energy as S OT (e„), where 

H ml {E) = -ie 2 J d^xd 3 y J j^^^^ m {^S F {x,m (2) 

and a self-mass counterterm needed to renormalize £ is understood to be included. If we 
set n = 4, carry out the d 3 k integration, and represent the Dirac-Coulomb propagator by a 
spectral decomposition, the above can be written 



sus) = <a i***v v e ' v ^';-" 1 m^m M«^Mff> . (3) 

J J 2n r \x — y\ E — k — e r (l — id) 



If we define 



,, e iVE 2 +i8\x-y\ 

g ijk i{E) = a d 3 x d 3 y — — — - — ^i{x)^ k {x) ipjiyh^iiv), ( 4 ) 

J \3C 1j \ 



\x-y 

then the self-energy can be compactly represented as 



i \ . 1" dko 9vmmv(ko) /r\ 

A; "V 2^^e v -k -e m (l-i5Y ^ 



It will be convenient below to also introduce the function 

g m {E) = a[d 3 x d 3 y sin (^-^l) ^(f) 7 ^ fc (f) faffi-ffrtf). (6) 
J \x y\ 

To carry out the numerical evaluation of the Lamb shift, a Wick rotation with k — > iu 
is performed. The resulting expression is purely real because imaginary parts present in the 
uj = — oo interval cancel against other imaginary parts in the u = -oo — interval. The 
imaginary part of the self-energy arises solely from the pole term, where a bound state pole 
in the first quadrant present when e v > e m is encircled during the Wick rotation. As we are 
interested in decays to the ground state, we will not consider imaginary parts of the energy 
arising from m being an excited state in Eq. (0). We introduce the convention that a refers 
to the ground state when there is no dependence on the magnetic quantum number (as is 
the case for the energy e a and the self-energy E aa (e a )), and b or c refers to the ground state 
when there is a dependence, with a sum over b or c running over the two possible values 
(±1/2) of the magnetic quantum number. We also define the lowest-order decay photon 
energy AE = e v — e a . It is important to emphasize that this energy differs from the actual 
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photon energy because the energy levels are shifted by radiative corrections: the effect of 
these shifts will be accounted for perturbatively below. The pole term is 



J pole 



X>»»(A£). (7) 



In calculations of the Lamb shift the real part of this is taken, but here we are concerned 
with the imaginary part, which gives the lowest-order decay rate 

T = -2j29vbbv{AE), (8) 

b 

which can be written CIS db partial wave expansion 

r = -87raA£ £ Jd z xj l {AEx)Y lm {n x )^ v {x^ ll ^ h {x) 

blm 

x Jd 3 y Jl (AEy)Y l * m (n y )My) 1 ^ v (y)- (9) 

Because Feynman gauge has been used for the self-energy calculation, this form of the decay 
rate is different from derivations which use the properties of the actual transverse photons 
that are emitted, but the result is the same because of gauge invariance. In Table [fl we 
present the lowest-order rates for the states 2si/2, 2px/ 2 , and 2p 3 / 2 to decay by one-photon 
emission to the ground state for Z = 5, 10, ... , 100. We now turn to the radiative corrections 
to these decay rates, which we define in terms of a function R(Za), 



r = r n 



l + -R(Za)]. (10) 
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Before presenting the exact calculation of R(Za), we calculate the leading contribution of 
order (Za) 2 for the p states using an effective field theory approach. We do not treat the 
more complicated s state correction, as the Ml decay is highly suppressed at low Z. 



III. PERTURB ATIVE APPROACH 



When Z is small, an expansion in Za converges rapidly. We present here the calculation 
of RiZa) to leading order (Za) 2 , which will serve as a check of the nonperturbative treat- 
ment presented in the next section. The radiative correction to the decay rate is obtained 
from the nonrelativistic form of quantum electrodynamics (QED) supplemented by one-loop 
corrections to electron form factors Fx, F 2 and the vacuum polarization. In the lowest order, 
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the decay rate of the 2P state in hydrogen-like atoms is 



-aE 3 
9 



(lS\f\2P) 



where E is the nonrelativistic limit of AE defined in the previous section, 

E = E { 2P)-E(lS) = ^l, 

with the nonrelativistic wave functions 

1 



>is 



Hp 



Ait 
1 1 



(mZa) 3/2 2e 



mZar 



(mZa) 



3/2 -mZar/2 



(mZar). 



(11) 
(12) 

(13) 
(14) 



4tt 2^6 

Natural units in which H — c — 1 are used here. Note that <p2p is normalized here in a 
nonstandard way, namely / d 3 r (f>2P • $2P — 1- Using the nonrelativistic matrix element 

1 256 1 



d = (lS\r\2P) 
Eq. (|TT|) gives the well-known decay rate 



x/6 81 mZa 



T = -ma(Za) 4 . 

The radiative corrections to this can be expressed as 

a „ 5r 5rf 
-i? = — =3 — + 2-. 
7r 1 o E a 



(15) 



(16) 



(17) 



When QED effects can be treated as local potentials, the calculation of radiative corrections 
is relatively simple. We illustrate this with the correction due to the presence of vacuum 
polarization, which is given in the nonrelativistic limit by a local interaction potential 

AZa 2 



SV 



15m 2 



(18) 



Corrections to the energy and wave function of the 2P state from SV do not contribute to 
R(Za) at the order of interest, but the potential shifts the IS energy by 

4ma(Za) 4 



6E(1S) 



15vr 



(19) 



which gives a contribution to R(Za) of 32/15 (Za) 2 . In addition the potential shifts the IS 
wave function by 

1 



(E - H)> 

%a{Za 2 ^ 



SV 



mZar 



15tt 



IS 



- — 1e — mZar H 

2 2mZar 



ln(2mZar) 



1 



{mZa) 3/2 , (20) 



which leads to a total contribution from vacuum polarization of 



39 

Ryp(Za) = (Za) 2 \ ±- + 
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{Zaf 
{Zaf 



Ami 

15 3 

41 



131 



90 



61 8 , , 
In - 

90 15 3 



0.524 347 



(21) 



We note the strong cancellation between the effect of the energy shift, which is automatically 
accounted for when experimental energies are used, and the perturbed orbital, which shows 
care needs to be taken when that approach is taken. We now turn to the more complex self- 
energy correction. The effect on the energy shift is well-known, coming from the self-energy 
part of the Lamb shift of the IS and 1P\ji states of 



6E(1S) 
6E(2P) 



ma 
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ma 

7T 



{Zaf 
(ZaY 



■lnfco(lS') 



1 

6 



lnfc (2P) 



where 



\nk (lS) = 2.984 128 556, 
ln£;o(2P) = -0.030 016 709. 



(22) 
(23) 

(24) 
(25) 



This energy shift contributes to the decay rate in accordance with Eq. f)17jl . However, the 
radiative corrections to the dipole matrix element are more difficult to obtain. We split 
this correction into three parts, 8d — bdh + 5dn + Sd>K, where 5di comes from low-energy 
photons, Sdn is the high-energy correction to the wave function, and Sdx is the correction 
to the dipole operator. Using nonrelativistic QED one derives the following expression for 
Sd L 

2a 



Sd L 



37rm 2 Jo 
Pi 



S 

+ (s 
s 



LUdiU$l[f(Lu)] 

1 



(26) 



T 3 



Pi 



H - E s + uo H — E P + u 

1 1 

- Pi — — — p t 



T ■ 

3 (H — E 



H - Ep + uo 
1 



Pi 



-lis 



H-E s + tu Pi {H-E s 
1 



Pi 



P, 



Pi 



Pi 



(H-E s + oo) 2 
1 

(H-Ep + lu) 2 



Pi 



S 



Pi 



(27) 



where e is assumed to be asymptotically large and i, j are vector coordinate indices. All 
matrix elements of f(uj) are calculated numerically using a finite difference representation of 
the nonrelativistic Hamiltonian. The integration with respect to u requires special treatment 
regarding linear and logarithmic in e terms. The large u> asymptotics is 



, A B 
f{u) « - + — + 



c 

5/2 



+ . . . 



(28) 



where 



.4 



B 



P 



S 



(h — E s y 

4nZa 5 3 (r) 



Pi) + [S 
1 

{h - E s y 



(H-Ep 



P; 



Pi 

a 

V 24 



-2d 



d (Zaf 4 21,i 



C = 2y/2B(Za). 



(29) 
(30) 
(31) 



The numerical integration in Eq. (J2T?|) . along a contour which omits poles from above or 
below, leads to the result 



2m 

Sd L = —{Ae + B\n 

377 



2e 



d(Zaf 17.759 359 



(32) 



.{Zaf. 

The term linear in e is dropped, and the logarithmic term is canceled by the contribution 
Sdn coming from large photon momenta. This latter contribution can be expressed in terms 
of an interaction potential 5V obtained from the one-loop electron form factors F\ and F 2 , 

Za 2 L-S 



SV = Za' 



2tt 



(33) 



It contributes to the energy shift in a way that has already been accounted for in Eqs. 
l23|) . but also gives corrections to the wave functions, and thus to the transition dipole 
moment 



8d 



H 



Za 2 
Za 2 

+ 



t -;-<*> 



S 



6 3 (r) 



(E s - H)> 



2tt 



S 



1 



L-S 



(E P - H)' r ? - 



P). 



(34) 



There is one more spin-dependent term which recently was discussed in Ref. It arises 
from the anomalous magnetic moment k correction to the dipole transition operator 



% u r — 



Am 



k 2 fx a. 



(35) 
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Its matrix element between S and P\ii states leads to a correction Sdx 

5d>K 



-d — 2 = -d^-(Za) 2 . 



Am " 32 vr v ' ^ 

With the help of Eq. (fTT|) . the sum ddL + Sdn + Sdx, together with energy shift contributions 
from Eqs. (J22| l23|) . leads finally to the result for the radiative correction to the decay rate 
of the 2Pi/2 state 



8,4 61 
-In 

3 3 18 



R 2 ^\Za) = {Za) 
A similar calculation for the 2P 3 / 2 state yields 



R 2 s T(Za) = (Zaf{ 



8,4 61 
-In 

3 3 18 



\n(Za)- 2 + 6.051 68 



\n(Za)- 2 + 5.984 36 



The coefficient of the logarithmic term is in agreement with 



121- 



(37) 



(38) 



IV. TWO-LOOP FORMALISM 

Following the approach given above to calculate radiative corrections to decay rates, we 
consider the imaginary part of the two- loop Lamb shift. We begin by considering the three 
self-energy diagrams of Fig. ^ leaving vacuum polarization for later. Expressions for the 
diagrams, which we refer to as overlap, nested, and reducible following the notation of Fox 
and Yennie were derived by Mills and Kroll |17| . and we now treat them in order. 

A. Overlap diagram 



The overlap diagram, Fig. [IJi, is given by 
S 4 ° = 



d n k d n l Q ik -(x-z) e il-(y-w) 

e 4 / d 3 x d 3 yd 3 z d 3 w / ■ 77 ^ - : - J9 : - i/) v (x) Y 



(27r) n (2ir) n k 2 + id P + id 
x S F (x, y; e v - k )YS F (y, z\ e v -k - k^Spiz, ™; e v - loh^ v (w). 



(39) 



As with the one-loop case, we introduce spectral representations for the electron propagators 
and carry out the d 3 k and d 3 l integrations to get 



^40 



^ r dk f dl 



Qvnmr (ko) Qmrnv 



2n (e v — ko — e m )(e v — k — l — e n )(e v — Iq — e r 



(40) 



where in this section we leave the factor (1 — iS) multiplying energies in the spectral rep- 
resentation of the electron propagator understood. We now consider Wick rotating both 
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ko — > ic^fc and Zo - > ^f- If no poles are passed, this again leads to a purely real expression. 
To get an imaginary part, at least one of the three denominators must involve encircling a 
pole. The middle denominator can have a pole when n = a and k + l = AE, but this 
corresponds to two-photon decay which we do not treat here. We then need consider only 
the two cases when either m = a and ko = AE or r = a and Zo — AE, which gives the 
expressions 

e iAE\x~z\ . ^n^ e ih{y-w) 



Ef = -4nia 2 Y, ! d 3 x d 3 y d 3 z d 3 w^— — — /- 

i J \X Z\ J I 



and 



\x-z\ J (2vr) n P + i5 
x i) v (x)^i) b {x) $b(y)'fS F (y, z\ e a - lo^Spfa w; e v - k)^ u ip v {w) (41) 



e iAE\y-w\ r d n k Q ik '(%-z) 



E 4 R ° = -Ania 2 y d 3 x d 3 y d 3 z d 3 w — — / 

b J \y-w\ J 



y — w\ J (2n) n k 2 + i5 
x ^ v (x)^S F {x, y; e v - k )YS F (y, z; e a - k )^^ b (z) ^(w^ip^w), (42) 

where we have "undone" the spectral representations of the electron propagator and kept 
either the d 3 k or d 3 l integration. 

We note at this point that these expressions are almost identical to expressions that arise 
in the treatment of screening corrections to the self-energy in lithiumlike ions (Eqs. (25, 27) 
in Ref. Q|), with the only difference being an overall minus sign and the fact that we are 
interested in the imaginary part here, while the real part was calculated in ^|. We were 
able then, with only slight modifications, to use code developed for the screening corrections 
in lithiumlike ions for the present calculation. Replacing 

e iAE\y-w\ _^ sin (A£|y- ^|) (43) 

and using the equality of the two terms gives the net result for the decay rate contribution 
from the overlap diagram we call IV, 

IV = -167ra 2 3V / d 3 xd 3 yd 3 zd 3 w sin ( AE \ S ~ ^ ^(f) 7 ^(f) 
b J \x-z\ 

[ d n k g**'^ - ^ 

x J j^yi k 2 + iS ^b{y)l u S F {y, z; e a - k^^S^z, w; e v - k G )^ v {w). (44) 
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B. Nested diagram 



The nested diagram, Fig. [TJd, is given by 

r r rl n b rl n l p ik-{x-w) il-(y-z) 

x S F (x, y; e v - k )YS F (y, z; e v - k - l Q )^ u S F {z, w; e v - k )^^ v (w). (45) 
This leads to the expression 

V"^ I I dlo 9vrmv{ko) 9mnnr(lo) (46) 

££ r J 2n J 2n (e v - k - e m )(e v - k - l - e n )(e v - k - e r ) 
We again consider Wick rotating both k — > iujk and Iq — > iuii, which gives a real result if 
no poles are passed. To get an imaginary part, at least one of the three denominators must 
encircles a pole, and once again, we omit poles arising from the middle denominator, which 
correspond to two-photon decay. We therefore need to consider only the ko Wick rotation, 
which has poles when ko = AE and either m = a or r = a. However, if both m and r are 
the ground state, a double pole is encountered. If the double pole is excluded, two terms 
result, 

S 4Ar = 5 f wfa(A£')Sfc T .(e a ) 

br €a ~ €r 

and 

Yj47V _ Qvbmv (AE) £ m fr (e a ) 

bm ta ~ €m 

These terms can be written in terms of perturbed orbitals. Specifically, if we define 

& (5) = £ W (5) 9vrbu ^ E) 5 mbmr , (49) 

fer ta Gr 

and 

^a(^) = £ i) m {z)^^2 <Wm> (50) 

fern Ca 6m 

then 

+ S 4 / = £ a5 (e a ) + E 5a (e a ). (51) 

Because the ground state self-energy is purely real, the only contribution to the decay rate 
comes from the imaginary part of these perturbed orbital terms. 
To treat the double pole, we set m = b and r = c, and find 

„47V • \ - f dk 

Qvcbv 

(k Q )H bc {e v - k Q ) 

hD 2vr (AE-ko + zSy • [bZ) 



be 
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Applying Cauchy's theorem and using the fact that the self-energy is diagonal in magnetic 
quantum numbers gives two derivative terms, 

b b 

C. One-Particle Reducible Diagram 

The final contribution to the two-loop self-energy comes from Fig. which breaks into 
two parts, a perturbed orbital term 

v-i \ ^vm (^l)) ^mv i^-v) /rA\ 

l po = Z ' ( 54 J 



and a derivative term 

dE vv (E) 



(55) 

E=e v 



dE 

The perturbed orbital term will have an imaginary part only if at least one pole term 
is present, as our analysis of the complex nature of S did not depend on the external 
wavef unctions, so long as they are real. This then leads to an imaginary contribution to the 
energy of 

Spo(a) = i[%vv(e v ) + S c „(e„)], (56) 

where 

U^ = ^ r {^ ~ 9vrbu{ i E \ (57) 

br tv €r 

and 

U% = E 9vbmv ^ ] . (58) 

bm 6v 6m 

The derivative term will lead to an imaginary part of the energy in two ways: in the 
first, we take the imaginary part of the first self-energy, which is of course associated with 
the lowest-order decay rate, and multiply it by the real part of the derivative of the valence 
self-energy. We combine this term with the first term of Eq. (J53j) to get the "derivative A" 
term, 

T dera = T [Ka(e a ) + &Kv{e v )}. (59) 

The second contribution is when the real part of the first self-energy multiplies the imaginary 
part of the derivative of the self-energy, which can be combined with the second part of 
Eq. (|53~|l to give 

^derv = To [^w( e u) ~ ^aa( e a)]- (60) 
11 



In our numerical analysis, we simply evaluate T' Q as one object. However, as can be seen 
by referring to Eq. Q, a multiplicative factor AE is present in the formula for IV If the 
derivative acts on this term, a contribution of T /AE would be present, as is the case in the 
formulas given by Barbieri and Sucher 

V. VACUUM POLARIZATION 

While the exact treatment of vacuum polarization is somewhat complicated, to order 
a(Za) 2 one needs to consider only the analog of the Is perturbed orbital. This is to be 
contrasted with the effective field theory discussion, in which both that perturbed orbital 
and an energy shift needed to be considered. While the effect of the energy shift is present 
in the exact calculation, which enters through the analog of the derivative A term, it is 
a peculiarity of the Feynman gauge that the \ow-Z behavior of r' is of order (Za) 4 : this 
arises through a cancellation between timelike and spacelike terms, which separately behave 
as (Za) 2 . Replacing the self-energy with the Uehling potential in the Is perturbed orbital 
gives numerical results that are consistent with Eq. ()21|) . 

VI. REARRANGEMENT FOR NUMERICAL EVALUATION 

In this section we perform further manipulations on the basic expressions for the two- 
loop self-energy that will allow an exact numerical evaluation. Beginning with the overlap 
term, we note that it is ultraviolet divergent. We deal with that divergence by considering 
Eq. (JUJ) with the bound state propagators replaced with free propagators, which leads to 
an expression we denote Tyi, 

rvi = - 2 -^%Jd*xJ dh sin( yjl~ ^ E MV-fMV J d 3 P 2 J d 3 Pl 

r d n k 1 1 1 

x ^'^S^ WTTs ^^l^Tn^A^f^/^ (61) 

where P2 = (e a , p<i) and p± = (e v , pi). If we define q = f>2 — Pi and 

J^Q) = J d'xMx^M^e- 1 ^, (62) 

this can be rewritten as 

12 



X ^ 2)7 V l-m ^ ^- l-m ^W- (63) 
Standard Feynman diagram techniques can now be used to write this as 



VI 



7T6 

2 



a 

+ 



2n 2 AE h 



£ fd 3 p 2 d 3 Pl J^q)5(\q\-AE) f P dp( l dxN^\n— 2 
£ J d 3 p 2 d 3 Pl J» b (q) 5(\q\ - AE) jf prfp ^ rfa; (64) 



where 



Att 2 AE , 



C = (47r) s / 2 r(l + e /2), 

Q = p[xpi + (l-x)p 2 ), 

A = px(m 2 — pi) + p(l — x){m 2 — pi) + Q 2 , 



and 



Nf, = Mfchvifa- Q + m)^i- Q + m)YMPi)- (65) 
We note that the momentum space form of the lowest-order decay rate is 

r = -^^E Jd 3 p 2 d 3 pMq)m\-^E)N ,. (66) 

The ultraviolet divergent term in Ty\ will be shown to cancel with derivative terms below, 
and the finite remainder is tabulated as Qy\ in the second columns of Tables HP- II VI where 
we adopt the convention 

T X = -T Q X . (67) 

7T 

We can now deal with an ultraviolet finite expression by evaluating IV — IYi. To numer- 
ically evaluate the subtracted form, we first carry out the Wick rotation ko — > ioj. If this 
passes no poles, it is straightforward to show that an expression we refer to as Ty 2 results, 
4a 2 , sm(AE\x - z\) r™, e^ 1 ^ 1 



T V2 = > / d A x d A y d A z d A w — ._ _ — - 3? duu 

ix ^ J \x — z\ Jo 



\y - w\ 



x %l) v (x)^^ a (x) ip a (y)-i u S F {y, z; e a - iu)^S F (z, w; e v - iu)Y^v(w), (68) 

where a subtraction of the same form with free electron propagators is understood. A kind 
of infrared divergence called a reference state singularity is present in the above, and is 
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regulated by taking e v — > e v (l — 5) and e a — > e a (l — 5), where 5 is chosen here to be 1CT 6 . 
Tv-2 has a logarithmic dependence on S which cancels with derivative terms and this is one of 
the checks used in the calculation. It is possible to combine the terms together to manifest 
the cancellation, but we have found it simpler to work with a small, finite value of 5, checking 
of course that the sum is unchanged when S is varied. We tabulate Qv2 in the third columns 
of Tables ITTl-llVl 

Finally, the Wick rotation picks up pole terms. To treat these, it is convenient to rewrite 
Eq. flSJ) as 

V -AO* \ ^ Qvmbn 

(AE) 

Qbnmv 

! J 2tt £^ (e a - k - e m )(e v - k - e n ) ' 

Because of the regularization procedure, the first term in the denominator has no poles, but 
the second does when e n < e v , which leads to the pole term 

(AE) 

i-QV3 = ^2^ T — Z — Z n ' ^ > 

bmn 6m ' Ev 6(1 6n 

where F n = 1 for positive energy states with e n < e v and F n = otherwise. The associated 
contribution Qvz is tabulated in the fourth columns of Tables ITT1 - IIVI 

Evaluation of the derivative A terms of Eq. (|59|) is similar to that of Qv2, as in both cases 
ultraviolet and reference state singularities are present. The same procedures are used to 
deal with this, namely a subtraction of a free-propagator term and use of the 5 regulator. 
The analog of Qvs is also present, although in this case it involves a double pole. Since we 
have discussed the evaluation of these derivative terms in some detail in a number of other 



papers (see, e.g., |1& l20||). here we simply combine the various finite effects into the terms 
Qsl2 and QsR2 in the tables, where Qsl2 refers to T>' vv and Qsr2 to S' aa . An ultraviolet 
divergent term in the free-propagator terms cancels the first term in Qyi- The perturbed 



orbital terms are evaluated using techniques for evaluation of the one-loop Lamb shift j21j . 
with Eq. tabulated as Qsli and Eq. (|5Tj) as Qsri- Finally, the derivative B term of 
Eq. (jnni) is tabulated as Qderb- 

VII. DISCUSSION 

The most numerically striking feature of the present calculation is the very large degree 
of cancellation present in the 2s Ml decays, which prohibits going below Z = 50. In the 
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lowest-order calculation, while using Feynman gauge gives the correct answer, a large cancel- 
lation between a timelike and spacelike contribution is present, leaving the highly suppressed 
(Za) 10 result shown in Table This cancellation is lost in individual contributions to the 
radiative correction, and is only restored in the sum. This strong cancellation in fact served 
as a useful test of the formulas and numerical methods. In the unlikely event that radiative 
corrections needed to be considered for Ml decays in hydrogenic ions with lower Z, the 
calculation would be better carried out in the Coulomb gauge. 

Turning to the 2p El decay rates, we note that, while less severe than for Ml decays, 
there is still considerable cancellation present between the various contributions, particularly 
at low Z. This is of course required by the fact that the Za expansion series for R(Za) has 
no constant term, but instead starts in order (Za) 2 . Again, the high degree of cancellation 
between contributing terms at low Z serves as a check of our numerical calculations, but 
in this case, we can also compare our low-Z results with the perturbation series. In Fig. |2l 
all-order results of Rse(Zoc) for 2pi/a and 2p 3 / 2 from Tables HTD and UVl are compared with 
analytic results from Eqs. (|H7|) and (|H8jl . It can be seen that all-order results do converge 
to the analytic results at low Z. In particular, the leading logarithmic term in the pertur- 
bation series, which is the same for both 2py 2 and 2p 3 / 2 , is good for Z = 1 and 2 only, 
while including the constant terms, which leads to the splitting of fine structure results, 
extends the validity of the perturbation series to Z = 7 or 8. This is typical of self-energy 
calculations where the Za series is known to converge very slowly except at very low Z and 
nonperturbative methods such as those shown here are needed for mid- to high-Z ions. 

In spite of the apparent agreement between the perturbation and all-order results shown 
in Fig. it should be noted that there are residual, unresolved discrepancies between them. 
By extending the accuracy of our all-order calculations for R$s(Za) to a level of ±0.0002 for 
low-Z ions, we were able to extract values of 6.67 and 6.62 for the constant terms of the 2pi/ 2 
and 2p 3 / 2 states, respectively. While these results are uncertain to ±0.20 due to the high 
degree of numerical cancellation at low Z when the leading logarithmic term is subtracted 
out, they are nevertheless different from the corresponding analytic values of 6.051 68 and 
5.984 36 from Eqs. (|37jl and (J38j) . Until this discrepancy is resolved, we would assign a 10% 
error to the constant terms, which should have negligible effect on R^(Za) anyway. 

While the decay rate corrections here are of intrinsic interest, the purpose of the present 
calculation is actually to serve as the first step in the calculation of PNC corrections. There 
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is interest in the parity forbidden transition 6S1/2 — > 7si/2 in cesium, which serves as one 
test of the electroweak part of the standard model of particle physics. A very large radia- 



tive correction has been found for this case 



22j . [23], but the only calculation using exact 



propagators that has so far been carried out was for the 2sx/2 — 2pi/2 matrix element for 



hydrogenlike ions 



2J]. This calculation had the advantage of being gauge invariant because 



of the degeneracy of the Dirac energies of the two states. The calculation carried out in 
this paper is also gauge invariant despite the differing energies of the n = 2 states and the 
ground state, and is generalizable to the PNC case. To carry out this generalization, the 
extra perturbation of the effect of Z boson exchange between the nucleus and electron must 
be added. 

An additional feature that must be dealt with for extending the present calculation to 
neutral cesium is the complication of dealing with a many-electron system. Because we use 
numerical Green's functions, there is no difficulty in using a modified Furry representation 
of QED in which the Coulomb potential is replaced with a model potential that incorporates 
the dominant effect of electron screening. However, the effect of the filled xenon-like core 
will have to be taken into account, which will lead to extra diagrams. We are at present 
setting up calculations of radiative corrections to allowed transitions in the alkalis, in which 
these issues will arise, with the next step being the inclusion of the effect of Z exchange. 

The challenges to experiment in testing the calculations presented here are considerable. 
The largest radiative corrections found here are those to Ml decays at high Z, with the 
correction at Z — 100 being 1.2%. For El decays, even the largest case, 2py 2 at Z = 100, 
has a radiative correction of only 0.2%. Rather than a direct measurement, experiments 
involving interference, such as the one discussed in Ref. [25], may be more promising. 

There is a radiative correction, even though very small at low Z, that is of particular 
interest. It involves the El decay of the 2pi/ 2 state in hydrogen. One approach to the 
determination of the Lamb shift as described in Ref. |l3j involves the measurement of the 
decay rate of the 2pi/ 2 state in hydrogen to very high accuracy. To interpret the experiment, 
Ref. Q used the following formula for the radiative correction 



R2 Pl/2 (Za)^ { y = T (Za) 



1 1 19 

- In(Za)- 2 - - \nk (2P) + \nk (lS) - - - - 



(71) 



which can be shown to be equivalent to the first term of Eq. ()17j) . However, as discussed 
in connection with the vacuum polarization contribution, using only the energy shift gives 
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answers in significant disagreement with using both parts of Eq. (JT7J). Our result, combining 
vacuum polarization with self-energy, is 



R 2pi JZa) 



8,4 61 
-In 

3 3 18. 



In(Za)- 2 + 6.576 03 j. (72) 



present work 

As noted earlier, we do have agreement with the logarithmic contribution found earlier by 



Ivanov and Karshenboim |12|. Using the hydrogenic value of To (including recoil through a 
factor of m r /m e ) 



r (2p 



1/2, 



Z=l 



6.264 942 3 x 10 s s" 1 , 



the radiatively corrected lifetimes of the 2pi/2 state for hydrogen are 

r(2pi /2 ) = 6.264 927 4 x 10 8 s^ 1 (present work), 
r(2pi /2 ) = 6.264 922 3 x 10 s s" 1 (Ref. |l2|), 
r(2p 1/2 ) = 6.264 881 2 x 10 8 s" 1 (Ref. |l3|). 



As indicated earlier, Ref. jl2j included only the logarithmic term, which was in agreement 
with our results, so the numerical difference shown above is due to the constant term in 
Eq. (J72j) . The difference is under 1 part per million (ppm), which corresponds to under 1 
kHz in the Lamb shift. However, there is a more significant 7 ppm difference with Ref. jl^ . 
which should play a significant role in the interpretation of that experiment. Of course, this 
is only relevant if ppm precision can be reached experimentally. Issues involved in reaching 
this extremely high accuracy, which we note is two orders of magnitude greater than found 

It n 

in positronium p, have been discussed by Hinds 26]. 
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(a) (b) (c) 

FIG. 1: Two-loop Lamb shift diagrams. 




5 10 15 20 25 

Z 

FIG. 2: Comparisons between all-order and perturbative results of R^{Za). Solid and dashed 
lines are 2px/ 2 an d 2p 3 / 2 results, respectively. Closed and open circles are 2p 1 / 2 and 2p 3 / 2 all-order 
results, respectively. 
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TABLE I: Lowest-order one-photon decay rates to the ground state for n = 2 states of hydrogenic 
ions in atomic units. Numbers in square brackets indicate powers of 10. The last column gives the 
nuclear fermi distribution parameter c in fermis. Conversion to the unit of s _1 is through 1 a.u. = 
4.134 137 x 10 16 s-\ 



z 


2si/ 2 


2pi/ 2 


2p 3 /2 


c 


5 


5.9038[-16] 


9.4779 [-6] 


9.4735 [-6] 


1.8104 


10 


6.0733[-13] 


1.5172[-4] 


1.5144[-4] 


2.9889 


15 


3.5262[-ll] 


7. 6868 [-4] 


7.6546 [-4] 


3.2752 


20 


6.3251[-10] 


2.4321[-3] 


2.4140[-3] 


3.7188 


25 


5.9680[-9] 


5.9461 [-3] 


5.8769[-3] 


4.0706 


30 


3.7551 [-8] 


1.2351[-2] 


1.2144[-2] 


4.4454 


40 


6.9521 [-7] 


3.9207[-2] 


3.8037[-2] 


4.9115 


50 


6.8431 [-6] 


9.6268[-2] 


8.8114[-2] 


5.4595 


60 


4.5463[-5] 


2.0100[-1] 


1.8737[-1] 


5.8270 


70 


2.3180[-4] 


3. 7535 [-1] 


3.4051[-1] 


6.2771 


80 


9.8091 [-4] 


6.4597[-l] 


5.6706[-l] 


6.6069 


90 


3.6293[-3] 


1.0440 


8.8110[-1] 


6.9264 


100 


1.2193[-2] 


1.6033 


1.2913 


7.1717 



TABLE II: Breakdown of contributions to R2 Sl/2 (Za). 



z 


Qvi 


Qv2 


QV3 


Qsli 


QsL2 


Qsm 


QSR2 


Qderb 


R(Za) 


50 


-8418.269 


-68461.495 


78211.835 


-10142.571 


11.673 


8593.159 


11.810 


192.382 


-1.476 


60 


-2446.426 


-14508.277 


17540.361 


-4213.698 


11.696 


3490.708 


11.893 


110.828 


-2.915 


70 


-824.257 


-4479.499 


5587.689 


-1990.665 


11.725 


1610.892 


11.995 


68.086 


-4.034 


80 


-305.499 


-1593.601 


2044.727 


-1034.946 


11.762 


817.408 


12.121 


43.683 


-4.345 


90 


-119.499 


-639.062 


834.541 


-580.832 


11.808 


446.644 


12.277 


28.790 


-5.333 


100 


-47.230 


-288.465 


373.852 


-347.601 


11.867 


259.536 


12.474 


19.246 


-6.321 
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TABLE III: Breakdown of contributions to R2 Pl/2 (Za). 



z 


Qvi 


Qv2 


QV3 


Qsli 


QSL2 


Qsm 


QSR2 


Qderb 


R(Za) 


5 


-9.551 


-36.706 


23.027 


-0.013 


11 


.626 


-0.021 


11.625 


0.000 


-0.014 


10 


-6.887 


-28.062 


11.756 


-0.044 


11 


.625 


-0.064 


11.630 


0.000 


-0.045 


15 


-5.403 


-25.850 


8.106 


-0.082 


11 


.627 


-0.119 


11.640 


-0.001 


-0.082 


20 


-4.409 


-25.047 


6.358 


-0.126 


11 


.631 


-0.183 


11.653 


-0.003 


-0.126 


25 


-3.685 


-24.732 


5.370 


-0.173 


11 


.635 


-0.252 


11.670 


-0.003 


-0.172 


30 


-3.131 


-24.632 


4.760 


-0.222 


11 


.641 


-0.325 


11.691 


-0.012 


-0.230 


40 


-2.344 


-24.662 


4.097 


-0.319 


11 


.655 


-0.473 


11.743 


-0.031 


-0.334 


50 


-1.822 


-24.806 


3.793 


-0.409 


11 


.673 


-0.623 


11.810 


-0.065 


-0.449 


60 


-1.464 


-24.984 


3.654 


-0.488 


11 


.697 


-0.770 


11.893 


-0.121 


-0.583 


70 


-1.219 


-25.166 


3.594 


-0.553 


11 


.727 


-0.909 


11.995 


-0.207 


-0.738 


80 


-1.057 


-25.349 


3.570 


-0.599 


11 


.764 


-1.037 


12.121 


-0.334 


-0.921 


90 


-0.962 


-25.533 


3.558 


-0.629 


11 


.812 


-1.146 


12.277 


-0.521 


-1.144 


100 


-0.928 


-25.726 


3.544 


-0.637 


11 


.876 


-1.229 


12.474 


-0.800 


-1.426 
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TABLE IV: Breakdown of contributions to R,2 P3/2 (Za). 



z 


Qvi 


Qv2 


QV3 


Qsli 


QsL2 


Qsri 


QSR2 


Qderb 


R(Za) 


5 


-9.551 


-36.484 


22.775 


0.001 


11 


,626 


-0.006 


11.625 


0.000 


-0.014 


10 


-6.888 


-27.653 


11.254 


0.000 


11 


,625 


-0.014 


11.630 


0.000 


-0.045 


15 


-5.406 


-25.318 


7.397 


-0.001 


11. 


.627 


-0.023 


11.640 


0.001 


-0.082 


20 


-4.413 


-24.381 


5.418 


-0.002 


11. 


.630 


-0.033 


11.653 


0.004 


-0.123 


25 


-3.690 


-23.946 


4.202 


-0.004 


11. 


,635 


-0.042 


11.670 


0.008 


-0.166 


30 


-3.138 


-23.732 


3.368 


-0.008 


11 


,640 


-0.050 


11.691 


0.014 


-0.215 


40 


-2.352 


-23.572 


2.278 


-0.022 


11. 


,656 


-0.063 


11.743 


0.037 


-0.295 


50 


-1.830 


-23.555 


1.579 


-0.047 


11. 


,675 


-0.075 


11.810 


0.082 


-0.361 


60 


-1.470 


-23.594 


1.082 


-0.087 


11. 


,698 


-0.087 


11.893 


0.144 


-0.421 


70 


-1.218 


-23.661 


0.710 


-0.149 


11. 


,726 


-0.102 


11.995 


0.242 


-0.457 


80 


-1.043 


-23.751 


0.424 


-0.237 


11. 


.758 


-0.119 


12.121 


0.383 


-0.464 


90 


-0.926 


-23.863 


0.204 


-0.360 


11. 


,794 


-0.140 


12.277 


0.585 


-0.429 


100 


-0.856 


-24.006 


0.042 


-0.527 


11. 


.833 


-0.160 


12.474 


0.877 


-0.323 
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